UNCLASSIFIED 


AD  400  284 


HepAoduced 
itf.  the 


ARMED  SERVICES  TECHNICAL  INFORMAnWI  AGENCY 
ARLINGTON  HAU  STATION 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 

Government  thereby  incurs  no  responsibility,  nor  any 
obligation  i&atsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  foimdated,  furnished,  or  in  any  way 
supplied  the  sedLd  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  ri^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  nay  in  any  way  be  related 
thereto. 


PAUL  WEIDLINGER,  Consulting  Engineer 

770  Lexington  Avenue,  New  York  21,  New  York 


PART 


FORCED  VIBRATIONS  OF  AN  ELASTIC  CIRCULAR 
CYLINDRICAL  BODY  OF  FINITE  LENGTH 
SUBMERGED  IN  AN  ACOUSTIC  FLUID 


-COMPUTATIONAL  PROCEDURES  AND  NUMERICAL  EXAMPLE 


by 


MELVIN  L.  BARON  and  ALVA  T.  MATTHEWS 


Office  of  Naval  Research 
Project  NR  064-464 
Contract  Nonr  3454(00)FBM 
Technical  Report  No.  2 


January  1963 


■1' 


FORCED  VIBRATIOH3  OF  AH  ELASTIC  CIRCULAR  CYLIWmiCAL  BOOT  OF  PIHITE  LEWOTB 

SUBMERGED  IW  AW  ACOUSTIC  FLUID. 


Part  II  -  Computational  Proe<tdur»»  and  Hunwrleal  Exaaipl<. 


TABLE  OF  CONTENTS. 


List  of  S]niibola . . . . 

Introduction  ..  . 

IX  CoD^tatlonal  Procedures,  nf(0... . . 

Ill  Con^utatlonal  Procedures  and  Nuoerleal  Exaiqple,  n  •  0  .  •  •  • 

a)  Computation  of  tiae  Real  and  Imaginary  Portions  of  the 
Coefficients  .  Poji  »  . 

\i)  Determination  of  the  Set  of  Simultaneous  Linear 

Equations  on  the  Source  Strengths  0^^^  ,  and  C^j^  • 

0 )  Solution  of  the  Set  of  Simultaneous  Linear  Equations 
on  the  Source  Strengths  0^^  ,  B^^  and  C^j^ . 

d)  Evaluation  of  the  Fluid  Pressure  Field  . 

XV  Discussion  and  Conclusions . . . 

i^endix  A  •  Evaluation  of  the  Real  and  XmaginarF  Portions  of  the 
Coefficients  and  their  Space 

Dsrivatives  .  . 

Appendix  B  -  Uagrasaatie  Summary  of  Coaputatlons . 


Page 

1 

5 

9 

17 

20 

29 

32 

39 

hi 


59 

67 


•I' 


r,«,* 


r,a,5 


U,Y,W 


* 

* 

C 


*1 

•*2 

0 

0,B,L 


1 

J 


‘1.8.3  . 
‘i.a.S 


L 

I 


■ 


LIST  or  SYMBOLS 


*) 


Cylindrical  eoordlnataa  locating  aourea  pointaf 
aea  Fig.  (A-1). 

Cylindrical  eoordlnataa  locating  field  point) 
see  Fig.  (A-1). 

Longitudinal)  tangential  and  radial  diaplaeeatenta 
of  the  cylindrical  ahell. 

Radius  of  cylinder. 

Radius  of  small  circle  approximation. 

Source  strength  coefficient  per  unit  of  circumferential 

p 

length  for  acoustic  fluid  (in  /sec.). 

Velocity  of  dilatational  waves  in  elastic  cylinder. 
Velocity  of  shear  waves  in  elastic  cylinder. 

Velocity  of  sound  in  water. 

Source  strength  coefficients  per  unit  of  cireun- 

P 

ferential  length  for  elastic  solid  (in  ). 

Source  band  index. 

Field  point  index. 

Coefficient  ^  evaluated  for  o  ■  Cj^  )  Cg  )  c  as 
indicated  in  text;  £  is  non-dimensional  form  ~  . 

Length  of  cylinder. 

Rumiber  of  banda  into  which  cylinder  is  divided. 

Rumiber  of  cireuaferential  waves  in  the  cylinder 
di^placeIBentB• 


*) 


Additional  syiiiols  are  defined  aa  they  occur  in  the  text 
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>  PresBure  in  fluid  (Ib/in^). 

•  Expansion  coefficient  of  externally  applied  normal 

p 

traction  to  cylinder  (ib/in  )• 

•  Dietence  between  field  point  and  source  pointy 

-  See  Ea.  (IXI-16). 

•  .  Distance  between  field  point  and  origin. 

•  Contribution  of  small  circle  of  distributed  sources^ 
see  Appendix  C  of  Reference  [1] . 

•>  Time. 

-  Coiqputatlonal  varleble. 

•  Integral  coefficients  of  source  strengths. 

•  Elasticity  constants  of  cylinder. 

-  Potential  function  associated  with  velocities  of 
acoustic  fluid. 

•  Potential  functions  associated  with  displacements  of 
elastic  cylinder. 

-  Angle  between  field  and  source  points,  see  Fig.  (A-l). 

>  Mass  density  of  fluid. 

•  Ron'dimensional  variables  in  radial  direction. 

•  Non-dimensional  variables  in  s  direction. 

•  Cylinder  stresses. 

•  Angle  locating  small  circle  of  distributed  sources, 
see  Fig.  (C-1). 
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*  Frequency  of  ribration. 

-  Slope  of  t  see  Wg.  (5)* 

r  or  •  appearing  aa  a  aupereerlpt  for  the  coefficients  »  ^nji 

denotes  differentiation  with  respect  to  the  particular  waridble  used. 

Dots  indicate  differentiation  with  respect  to  tiise. 
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1  IMTRODUCTIOH 


The  theoretical  work  and  analyticeil  procedures  leading  to  the  evaluation 
of  the  pressure  and  velocity  fields  produced  In  an  acoustic  fluid  by  the  forced 
vibrations  of  an  elastic  circular  cylindrical  body  of  finite  length  vhlch  is 
submerged  In  the  fluid  have  been  presented  in  Reference  [1]  The  application 
of  the  theory  to  problems  of  practical  interest  Is  discussed  In  the  present 
report  and  numerical  results  are  presented  for  a  sample  axl-symmetrlcal  problem, 
l.ei  a  ease  In  vhlch  the  cylinder  excitation  is  Independent  of  9  and  the  response 
is  given  by  the  n  ■  0  component  only.  For  convenience,  specific  formulas  In 
computational  form  are  given  and  the  major  confutations  are  summarized  by  flow 
charts.  In  the  sections  vhlch  follow,  it  Is  assumed  that  the  reader  Is  familiar 
with  the  material  contained  In  [1]  and  that  all  unpreflxed  formiU.a  numbers 
and  section  designations  refer  to  the  material  In  that  report. 

A  potential  theory  approach  was  used  In  [1]  to  evaluate  the  pressure  and 
velocity  fields  In  the  fluid  due  to  the  tlme>harmonic  excitation  of  a  submerged 
cylindrical  body  of  finite  length.  The  stresses  and  velocities  In  the  elastic 
cylinder  vere  expressed  In  terms  of  three  displacement  potential  functions,  each 
of  which  satisfies  the  wave  equation.  Similarly,  the  corresponding  fluid  quantl* 
ties  vere  expressed  in  terms  of  a  single  fluid  velocity  potential.  Each  of  the 
four  potential  functions  vere  considered  to  be  caused  by  a  group  of  slnf le  sources 
of  unknown  strength  which  were  distributed  over  the  boundaries  of  the  elastic  body 
and  the  fluid  surface  at  the  eyllnder>fluld  interface.  A  finite  difference  approach 

"Forced  Vibrations  of  an  Elastic  Circular  Cylindrical  Body  of  Finite  Length 
Submerged  in  an  Acoustic  Fluid”,  by  M.L.  Baron,  A.T.  Matthews  and  B.H.  Bleleh, 
Paul  Veidlinger,  Consulting  Engineer,  Office  of  Naval  Research  Project 
]|R064>^A,  Contract  3^$H00)FBM,  Technical  Report  No.  1,  Jline  I962. 

2) 

'  The  numbering  of  fomulas  in  this  report  will  be  prefixed  by  the  aeotion  inad>er 
a.g.  Iq. 
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vu  used  In  vhlch  the  boundary  of  the  fluid- cylinder  Interface  vas  divided  Into 
a  series  of  bands.  On  each  band,  the  unknovn  source  strengths  were  considered  to 
be  constant.  Conditions  on  the  stresses  and  velocities  at  tho  fluid-elastic  body 
Interface  gave  rise  to  a  system  of  simultaneous  linear  algebraic  equations  on  the 
source  strengths.  The  coefficients  In  these  equations  are  definite  Integrals 
vhlch  are  evaluated  nximerlcally  by  suitable  quadrature  formulas  for  a  given 
geometry  and  forcing  frequency. 

In  order  to  solve  the  coupled  forced  vibrations  problem  for  arbitrary 
excitations,  the  tlme-harmonlc  exciting  forces  were  expanded  Into  a  Fourier  series 
In  6  around  the  circumference  of  the  cylinder.  In  this  manner,  each  term  corres¬ 
ponding  to  n,  the  number  of  circumferential  waves  In  the  particular  coirq>onent, 
could  be  treated  separately. 

The  solution  of  these  problems  Inherently  Involves  a  major  computational 
effort  directed  towards  the  evaluation  of  large  systems  of  linear  simultaneous 
algebraic  equations  on  the  source  strength  coefficients.  If  the  cylinder-fluid 
Interface  Is  divided  Into  N  bands,  sets  of  8n  (n  ^  O)  and  6n  (n  ■  0}  simultaneous 
equations  are  obtained.  It  will  be  shown  that  In  each  case,  the  computations  can 
be  reduced  to  the  evaluation  of  a  nuinber  of  systems  of  2N  equations.  For  the 
application  of  the  theory  to  practical  problems  In  which  as  mai^y  sis  50  or  mors 
bands  might  be  considered,  large  electronic  cosputers  capable  of  solving  systems 
of  one  to  two  hundred  linear  equations  would  be  required. 

Section  XZ  of  this  report  gives  the  procedures  and  detailed  formulas  for  the 
evaluation  of  the  n  ^  0  response  exponents.  Section  III  gives  similar  information 
for  the  axl-symmetrical  ease,  n  ■  0,  and  Includes  sn  Illustrative  numerical  exeople 
for  a  ease  with  the  parameter  ratios  ^  ■  2  and  ~  > 


2.01 


Xt  it  ftlt  that  the  present  theory  may  find  useful  applications  in  evaluating 


the  response  in  a  fluid  due  to  the  harmonic  excitation  of  large  transducers t 
additioni  current  work  on  the  extension  of  this  approach  to  the  case  of  thick 
vaUed  elastic  shells  submerged  in  an  acoustic  fluid,  is  under  way. 


In 
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II  COMPUTATIONAL  PROCEDURES,  n  i<  0. 

Tbit  leetlon  describes  the  computational  procedures  for  the  evaluation  of 
tbe  source  strength  coefficients  on  tbe  cylinder-fluid  interface  and  tbe  sub¬ 
sequent  evaluation  of  tbe  pressure  and  velocity  fields  In  tbe  acoustic  fluid. 

The  formulas  on  the  real  and  imaginary  parts  of  the  source  strength  coefficients 
are  presented  In  a  form  particularly  suitable  for  computations.  Tbe  cotiq>utatlonal 
procedure  Is  conveniently  divided  Into  three  major  portions  s  a)  Evaluation 
of  the  coefficients  and  and  their  space  derivatives,  using 

Eqs.  (29),  (30)  and  Eqs.  (U^)-(32);  b)  Solution  of  Eqs.  (80)-(67)  for  the 
source  strength  coefficients  and  C^^j^  j  and  c)  Evaluation  of 

the  pressure  and/or  the  velocity  field  In  the  fluid,  using  Eqs.  (67)«(69)  and 
Eqs.  (89)- (89a). 

For  a  particular  set  of  input  parameters,  l.e.  the  cylinder  geometryj  the 
elastic  constants  of  the  cylinder  material;  the  fluid  constants;  and  tbe  space 
distribution  of  the  tlme-barmonlc  excitation,  the  cylindrical  shell  Is  divided 
Into  a  number  N  of  subdivisions  and  the  computations  are  started.  The  number  N 
of  bands  on  the  shell-fluid  Interface  whlcn  are  chosen  for  any  particular  problem 
vlU  depend  essentially  on  the  complexity  of  the  applied  excitation  and  the  accuTM] 
required  In  the  results. 


The  evaluation  of  the  coefficients  ,  y^j^  and  their  space  deri¬ 

vatives  have  been  described  In  detail  for  the  n  ■  0  ease  In  Section  (IV)  of 
Reference  [1]  and  no  further  amplification  Is  required  here,  nie  conputatlonal 
formulas  for  the  real  and  Imaginary  parts  of  these  coefficients  for  n  ft  0  are 
given  In  ^endix  A  of  this  paper  for  the  tvo  cases,  1  ■  J  and  1  y*  J. 


3) 


See  Reference  [1],  Section  IV,  Pg.  37  ff 
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The  second  major  computational  effort  is  the  evaluation  of  the  complex 

source  strengths  0^^  ,  f  and  .  With  four  such  complex  unknown 

coefficients  for  each  hand  N  on  the  cylinder-fluid  Interfacei  this  leads  to  a 

system  of  8n  simultaneous  linear  algebraic  equations  on  the  real  and  imaginary 

4) 

parts  of  the  source  strengths  For  the  axleymmetrleel  case  n  «  0,  all  the 
coefficients  are  zero  and  a  system  of  6w  equations  is  obtained. 


The  source  strengths  and  the  coefficients  and  their 
space  derivatives  appearing  In  Eqs.  (8o}-(87)  are  written  In  terms  of  their 
real  end  Imaginary  parts,  e.g. 


®nl  "  ®nl  ^  °nl 


“nl  "  ®nl  *  ^  “nl 


and  the  real  and  Imaginary  parts  of  each  of  Eqs.  (80)-(87)  are  equated,  thus 

5) 

leading  to  the  following  set  of  equations  in  computational  form 

Point  J  on  the  surface  r  «  a 
R 

"•  V  -IK*!  -  *  *.^53  •  ‘»1*3  *  'nl^  * 

5.154} 


1-1 


+  c_ 


-  p 


rr,nj 


(II-l) 


r-a 


4) 

'  Recent  developments  In  the  programming  of  systems  of  simultaneous  linear 
algebraic  equations  with  con^lex  unknowns  may  allow  the  direct  evaluation 
of  the  eoiiq>lex  source  strengths  and  without  breaking  them 

Into  their  teal  and  Imaginary  parts.  !IMs  would  reduce  the  problem  to  the 
solution  of  systems  of  4N  (n  M  0}  and  3N  (n  ■  0)  equations  respectively. 

In  the  formulas  which  follow,  the  subscript  njl  in  each  of  the  d  coefficients 
will  be  understood,  e.g,  dg  ■ 
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n 

^  ‘Z{®nl*l  *  *  ®nl‘®2  ^  **111^3  *  ^nl*3  ’ 

5„i5u} 


*  ■  0 
r-a 


H 

*'*  ®r*  “Z{®»i^5  ■  °nl^5  *  ®nl‘^6  '  ®nl‘*6  *  “  ^niV  *  ° 

1-1 

N 

IB  c„  -  ♦  fi^adg  +  l^alg  4  L^d^  4  -  0 

1-1  r-a 

N 

®r0  “  *  ®nl®^9  "  ^nl^^9  *  ^nl^lO  "  ^ni^o}  “  ° 

1-1  r-a 

N 

^  "rt  ■  IK’s  *  *nl’8  *  ♦  ^nl\o  *  ^nl^oj  ’  » 

1-1  r-a 

Re  *  -^{-  4  a^ap”j  4  H^ea”^  -  f  - 

1-1 


a  c 

=^0  5  ^  nl  -r  nl  "r  1  ^ 

^ni  r  ^nji  *  to  ^ajl  “  to  ^ajlj 

r-a 


'■  ‘  -  -  5.a»Sjl  •  ^<1  ‘  ®ni*5Sj‘l  *  ?  5.J1 

■ 


C  c 

y  n  ■  _nl  ye  j^ 

“  “nl  r  ^njl  *  to  ”njl  ’’  to 


(II-2) 

(11-3) 

(II-*») 

(11-5) 

(11-6) 

(11-7) 


(11-8) 
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Point  J  on  end  surface b,  z  ■  0  and  g  ■  L 
N 


Raff  ,  -  O.d, ,  +  fi  .a3. .  -  fi  .ad,  _  +  C.d,  +C.3,  V  ■?  , 

SB  Zj  L  ^  ^  ^  ^  ni  4.J  BB,nJ 

!■!  zvO|L 


(n-9) 


N 

'bz  "  ■"  "■  “nl^^l2  “  Ki.\  *  °nl*4  }  "  ° 

!■!  ZrnOmlt 


(II-IO) 


Ra  ff  •  0  Apply  Eq«  (11-3)  at  point, 
z** 

In  ff  -  0  Apply  Eq,.  (II-U)  at  point. 

*  * 


N 

«•  »rf  ■  Zi’nl’ls  -  Vl3  *  -  "“nrtk  -  *  V^l}  ’  » 

1-1  b-0,L 


(ZI-U) 


n 

IK'13  *  '“nl^l3  *  *  “al'^1'*  '  ’  ‘nl5n;i}  ’  » 

1-1  z-O.L 


(11-12) 


^  "Z{’  “  ®nl®ajl  *  Ki^^*  u>  ^njl  ‘  id  ^njl}  "  ° 

1-1  z-0,L 


(11-13) 


n 

I"  ‘  'HvSji  -  -  8.i*3i5  * 

1-1 


e  c 

nl  “z  nl  -Z 
15  ■*■  CD  ’’njl  ■*■  CD  ^nji 


} 


0  . 


b-O^L 


(11-14) 
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Expressions  for  the  coefflclonte  3^  ond  are  given  in  Table  II-l  and 
the  oxpresoions  for  the  real  and  imaginary  porta  of  the  coefficients  t 

^nU  ^  ^nji  their  space  derivatives  are  given  in  Appendix  A  of  this  report 
[Bee  Eq.  (A-6)  -  (A-33)]. 


The  application  of  Eqa.  (Il-l)-(II-lJt)  at  the  points  on  the  fluid- cylinder 
Interface  yields  the  set  of  3N  simultaneous  linear  equations  on  the  soxirce  strengths. 
These  equations  con  be  solved  by  either  of  two  approaches: 

1)  a  direct  solution  of  the  system  of  6n  equations  on  a  large  electronic  computer,  o 

2)  a  reduction  of  the  system  of  8n  equations  to  a  series  of  2N  systems  end  their 
subsequent  solution  on  on  electronic  computer.  These  alternatives  will  be 
Illustrated  In  detail  In  Section  III  of  this  report  In  which  the  illustrative 
problem  for  an  axleymmetrlcal  problem,  n  ■  0,  Is  presented. 

Once  the  source  strength  coefficients  are  evaluated,  the  complex  fluid  pressure 
at  a  field  point  J  is  obtained  by  direct  substitution  of  the  source  strength  coef¬ 
ficients  Into  Eqs.  (89)“ (890'): 

N 

l<ut 


r*  y 

K' 


^ni^njl 


^nl’^nji 


) 


+  1 


(^nl^njl 


'  ^nl’' 


V"1 

M 


(11-15) 


The  far  field  pressinres  In  the  fluid  can  be  evaluated  from  the  simplified 
asymptotic  expression,  Eq.  (89a): 


lo)  R 


..  -  [-  - 


p.  cos  o.  ^  lk,5  sin  C 
- 2 -  )  • 


(kjr  cos  E)]  t— 


(11-16) 


-lu 


The  uao  of  this  asymptotic  expression  cT«BLtIy  simplifies  the  numerical 

computation  of  the  far  field  fluid  pressure,  since  it  eliminates  the  computation 

of  the  Integral  coefficients  for  each  point  along  a  specific  ray  In  the  fluid. 

An  additional  convenience  Is  tbs  fact  that  once  the  pressure  p_  .  has  been  evaluated 

n,  j 

from  Eq.  (Il>l6)  for  a  given  point  on  the  ray  specified  by  and  it  can 

easily  be  evaluated  for  any  other  point  on  that  ray  by  changing  the  scale  factor 
io)  R. 


,  [See  Figure  ZII»3l* 
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Table  Z  •  Evaluation  of  the  Coefficients  and 


6) 


; 

S  ■  [f  ■  7  *»ji] 

; 

*1.  ■  V  'ill 

^5  ■  ^n!l 

• 

j 

j  2n  -  2n  -r 

*8  •  7  *nji  *  T  ®njl 

j  2n  .  2n  rfa 

S  ’  y2  *njl  *  t  ^nji 

} 

2 

2 

; 

; 

2 

»  -  ff'*  *  22L  a 

T.4  r  ^njl  *  2  ^njl 

6) 


The  fomulas  for  the  coefficients  dj^  are  obtained  by  replacing  the 
single  bar  yalues  by  the  corresponding  double  bar  values  In  each  case 
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m  COMPUTATIOWAL  PROCEDURES  AND  NUICTICAL  EXA>gLE,  n  -  0. 


For  illustrative  pvurposes,  the  following  axl symmetrical  problem  Is  presented. 
An  elastic  cylinder  of  radius  a  and  length  L,  Immersed  In  an  Infinite  acoustic 
fluid,  undergoes  an  electromagnetlcally  Induced  time  •  harmonic  uniform  strain, 
c  ■  e  0^*°^  •  In  the  axial  direction,  while  the  radial  and  circumferential 
strains  and  Sgg  are  kept  equal  to  zero  (Fig.  Ill -la).  By  superposition,  the 
pressure  field  that  la  produced  In  the  fluid  by  the  straining  of  the  cylinder 
will  be  eqxilvalent  to  the  fluid  pressures  produced  by  a  set  of  fictitious  surface 
tractions  which  sure  applied  to  the  solid  cylinder  In  the  ratio 

o_ 


rr 


zz 


X.  +  2n  1  -  V 


(a) 


To  Illustrate  this,  consider  first,  the  cylinder  \mder  the  action  of  applied 
surface  tractions  and  o..  (Fig.  Ill-lb).  These  tractions  ere  chosen  so  as 
to  bring  the  cylinder  back  to  Its  original  unstrained  state,  l.e. 


rr  o 


Icut 


♦  2ti)  s^e 


luft 


(b) 


Finally,  a  set  of  surface  tractions  which  are  equal  and  opposite  to  those  of 
Bq.  (b)  are  applied  to  the  cylinder  (Fig.  XII-lc). 


The  Bxq^erpoBltlon  of  the  three  states  of  stress  of  Fig.  ( XZI-1)  shows  that 

the  pressure  field  which  is  produced  In  the  fluid  by  the  uniform  straining  of 

the  cylinder,  e  •  equivalent  to  that  which  Is  produced  by  the 

zz  o 

fictitious  surface  tractions  of  Fig.  (ZlZ-lc),  namely 


loot 


(\  +  2ji)  s^jC 


ICDt 


(c) 


-10- 


For  the  llluBtrntive  problem,  the  pariii,;:tcro  v  ■  and  ■  10^  have 

been  chooon.  The  problem  thuo  conslato  of  the  evaluation  of  the  prosoure  field 

produced  In  the  fluid  by  the  fictitious  time  harmonic  tractlona  e  ■  10^  e^^ 

rr 

and  c  •  3(10^)  e^^  which  ore  applied  on  the  boundary  surfaces  of  the  elastic 

%% 

cylinder. 


The  following  parameters  are  uned  In  the  numerical  example: 


Elastic  Cylinder 

2.  Acoustic  Fluid 

L  ■  2a 

w^  -  62.5  Ib/ft^ 

V  •  ^  (l.e.  K  * 

Cj  ■  5000  ft/cQC 

H  -  12(10^)  Ib/in^ 
w  «  O.2O33  Ib/in^ 

3.  Corputatlonol  Geometry  and  Pressure  Loading  [See  Fig. 

[6  side  bands  on  the  surface  r  ■  a 

6  ring  bands  on  the  surface  s  ■  0 

6  ring  bands  on  the  eurf ace  z  ■  L] . 


N  -  l8  bonds. 


k.  Symmetrical  Loading:  ■  (10^)  6^“^^ 


J  on  the  surface  r  >  a 


P  •  3(10^)  J  on  surfaces  z  •  0  and  z  ■  L 

zz,oj 


With  the  input  values  chosen,  the  propagation  velocities  of  the  pressure 
and  shear  waves  in  the  medium  become  respectively 


h  "  ^  “  18,W0  ft/sec. 

Cg  ■  ^  ^  ■  10,660  ft/sec. 


Cj  (fluid)  -  5,000  ft/sec. 


r- 

) 
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Tb«  following  values  of  the  noudimenslonal  parameters  R  m  were  used  In  the 

“  ®n 

eoiiq;>utatlonsi 

-  0.5444 

-  0.9425 

ic^  -  2.0100 

For  exaaqple,  a  possible  combination  of  forcing  frequency  to  and  cylinder  radius 
for  which  the  numerical  results  would  apply  Is  to  ■  400  cps  and  a  ■  4  ft. 

The  cylinder  geometry  for  the  choice  of  R  >•  I8  bands  Is  shown  In  Fig.  (ZZI>2). 
Zt  should  be  enphaslsed  that  this  value  of  N  was  chosen  only  as  a  reasonable  com¬ 
putational  gaometry  to  Illustrate  the  general  method  in  a  sensible  manner.  Zt 
f  does  not  necessarily  represent,  however,  the  number  of  bands  which  might  be 

required  to  obtain  physically  meaningful  results  of  sufficient  accurMy  for 
plication  to  practical  problems.  Zn  the  latter  case,  a  larger  number  of 
bands  R  would  generally  be  required. 


J 


•20' 


b)  Computation  of  the  Real  and  Imaginary  Portions  of  the  Coefflclente 
*oJl  *  ^ojl  *  ^ojl  their  Space  Derivatives. 


Some  representative  numerical  values  of  the  coefficients  a  ..  .  a  ..  .  their 

OJl  Ojl 

space  derivatives,  and  the  corresponding  (quantities  for  and  are  given 
in  this  section  for  the  Illustrative  problem.  Recalling  that  each  coefficient 
Is  vrltten  In  terms  of  a  real  and  Imaginary  part  e.g. 

“ojl  -  »0J1  ♦  ^  “ojl  ' 

analytical  e3g;>ressions  for  these  (quantities  are  given  In  Appendix  B  and  Appendix  C 
of  Reference  [!}. 


For  the  cases  In  which  1  »•  J,  the  quantities  *nd  their 

space  derivatives  are  evaluated  by  numerical  Integration  [Simpson's  Rule],  using 


Eqs.  (B-7)«(B-1$)  of  Reference  [1].  Sharp  peaks  which  occur  in  the  Integrands 
for  values  of  t  <  necessitated  the  use  of  a  variable  spacing  In  the  numerical 


Integration  as  follows: 


t(degreet) 

^t( degrees) 

0  -  8 

1 

8  •  20 

3 

20  •  1»0 

5 

Uo  •  6o 

10 

60  •  180 

20 

A  similar  numerical  Integration  was  used  to  evaluate  the  quantities  , 
^ojl  *  ^oji  derivatives,  Bqs.  (B>l6)-(B->24)  of  Reference  [1].  For 

each  quantity,  an  integration  spacing  of  ^  ■  16  degrees  was  used. 


21' 


Tor  ■mall  values  of  the  variable  Xi  eoiqputatlooal  accuracy  required  that  eertala 
of  the  trigonometric  eiqtreaalona  appearing  In  tbeae  equations  he  expanded  Into 
paver  series  so  that  the  leading  terms  in  the  Integrand  could  he  cancelled 
axialytlcally.  Expressions  for  the  evaluation  of  these  coefficients  are  obtained 
hy  setting  n  ■  0  In  Eqs>  (A-27)*(A>3?)  of  this  paper.  Equations  (A-36)'‘(A*>38) 
give  the  analytical  expressions  and  the  range  of  application  for  the  functions 
appearing  In  the  integrands.  It  should  he  noted  that  these  ranges  vlU  vary 
for  different  problems  and  that  they  must  he  determined  separately  In  each  case 
hy  a  series  of  trial  Integrations  using  the  expressions  given  In  Eqs.  (A-36)-(A'>38). 

For  the  case  1  «  J  in  which  the  field  points  and  the  source  points  are  on 
the  same  hand  special  formulas  are  required  for  the  evaluation  of  the  coefficients 
*oJl  *  ^ojl  '  ^ojl  hhelr  derivatives  because  of  an  infinite  discontinuity  In 
their  Integrands  at  the  value  defined  hy  R  ■  0.  To  evaluate  these  Improper  (hut 
convergent)  Integrals,  Eqs.  (126)-(127)  and  Eqs.  (C»4)-(C-8)  of  Reference  [1]  are 
used.  The  values  of  the  limit  $  In  Eqs.  (126)- (127)  are  evaluated  separately  hy 
trial  Integrations  in  eauih  case. 

Table  ZZ  gives  the  nondlmenslonal  values  of  the  quantities  p,  p  and  (t-() 
for  the  calculation  of  the  real  and  Imaginary  parts  of  the  coefficients  , 

and  and  their  derivatives  In  the  illustrative  problem.  While  It  Is 
not  practical  to  list  the  values  of  these  coefficients  for  all  combinations  of 
1  and  J,  some  representative  values  of  these  quantities  for  the  two  cases  J  ■  1, 

I  m  6  and  J  •  i  ■  9  *3^*  given  In  Table  III  for  the  convenience  of  those  readers 
vho  Bwy  vlsh  to  vork  through  a  saovle  computation.  Representative  values  for  the 
evaluation  of  the  real  portions  of  the  coefficients  and  Ita  pertinent  deriva¬ 
tives  are  given  in  Table  ZV  for  the  two  eases,  i  ■  J  •  5  (top  hand)  and  i  ■  J  «  7 
(side  hand)* 
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Table  III  -  Representative  Values  of  ,  7^,^  and  their  Space 


Derivatives  - 


CASE  X 


rtl  z 
a. 


Real  Portion 
(single  bar) 

Imaginary  Portion  1 
(double  bar) 

-.4858500 

+.2274170 

-.7521002 

+.OO38O7410 

+.6153040 

-.02273441 

+  3.910194 

-.0002299027 

+  47.68267 

+.001377773 

+  4.430315 

-.0002201813 

CASE  II 


J  -  3,  1  -  9 


Real  Portion  Imaginary  Portion 

(single  bar)  (double  bar) 


♦  2.07U596  -.02138003 


-.2860010 


•f. 2421691 


+.03017070 


+.04672338 


-.2926517 


-.O9OII3O7 


-.2010604 


+.2478359 


-.02119497 


-.01027833 


+.0005205666 


-.02437010 


+.003792008 


+.3522680 


-.09819865 


-.04451797 


+.007139976 


+.007779862 


-.1025629 


+.05333224 


-.3542777 

+.3262738 

-.7726809 

+.01807612 

+.6556329 

-.09997833 

+  3.9»*1736 

-.003183777 

♦  47.68707 

+.01903262 

♦  4.533949 

-.002781721 

+  2.119123 


■.08195736 


-.1225355 


+.2918503 


+.03667484 


+.03799412 


-.3136065 


-.09717981 


-.3392487 


m 

-,06006077 

+.2492022 

+.3OU386 

+.1673313 

-.8155342 

+.1005983 

+.1056342 

-.4913962 

m 

+.4507052 

-.3588395 

-.1396592 

-.1167767 

■87 


Tabl*  ZV  •  Reprts^nttttlv*  Valuet  of  and  ita  Spaet  Darlvativaa  — >  1  ■  j 


1  -  J  -  5 

1  -  J  -  7 

• 

2  -  -uflo 

m 

J  -  .1750 

Xa  ■  6.00 

Xa  -  2.97 

m 

a 

..6800 

-.63758 

5. 

-.3540 

-.2590 

aS* 

+  3*00 

0 

aS* 

a 

+  3*00 

0 

aSF 

+.34636 

+  1.8215 

as' 

a 

0 

+  1.4850 

.v* 

0 

0 

•2§” 

a 

0 

0 

.V“ 

+  59*820 

♦  U.401 

a3gr*« 

a 

0 

0 

a2a” 

•  20.582 

Or 

-  25.4o4 

mmm 

•V" 

-.889U 

... 

»3g*** 

a 

-.689U 

... 

.3s“*- 

•  •• 

0 

a 

mmm 

0 

... 

-  5*4126 

,2^ 

m 

— 

•  8.50 
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b)  Determination  of  the  8«t  of  Simultaneous  Linear  Equations  on  th»  Souro 
Strengtha  0^^^  ,  and  . 


As  for  the  caae  n  >«  0«  the  aourea  atrengtha  and  the  eoefflelenta  , 

7  ..  and  their  apace  derlvatlvea  are  written  In  terma  of  their  real  and  Imaginary 
ojl 

parta,  e4g. 


01 


'ol 


"ol 


"ol 


■»>  1  a 


ol 


and  the  reel  and  Imaginary  parta  of  Bqa.  (T17)>(122}  are  equated,  thua  leaidlng 
to  the  following  aet  of  equatlona  In  computational  form: 


Point  J  on  the  aurface  r  ■  a 
R 


»•  ‘rr  ■  *  °»A  *  'ol*^  •  *  'oA  ’  'oA}  ’  'rr,.J 


1-1 

H 


r-a 


•“  •«  ■  £{'oA  ‘  A  •  hA  *  'oA}  •  0 


1-1 

R 


r-a 


'r.  “  I{®ol^5  “  °ol*5  ®ol^^6  '  “ol^s}  ’  ® 


(III-3 


1-1 

R 


r-a 


^  •«  "  Z{®ol*5  ®ol®5  ®ol*^6  “ol'^^e}  •  ® 


(III-4 


1-1 

R 


r-a 


1-1 


r-a 


(111-5 


30- 


n 

1-1 


6  c 

rt  .  ol  rr  .  ol  j 


+  -7f  7;w  +  ~f  ^,4}  -  0  . 


ojl  (0  'ojl  0)  ojlj 


r-a 


(III-6) 


Point  J  on  end  eurfaceB,  x  -  0  and  t  ■  L- 
N 


*.*  -I{®oi^u  -  °oAi  ♦  ^oi^^ia  “  ^oi'^^ia  ♦  hiK  *  ^oi^u} 
1-1  t 


z-0|L 


iz,oi 


(111.7) 


n 

^  S.  ■  Z{*oAl  *  »ol*U  ‘  ‘  ■  0  »“-»> 


1-3 

Re 

s  -  0 
r* 

Apply 

Eq 

la 

9  -0 

ra 

Apply 

Eq, 

N 

Be 

•■a- 

A  ■* 
®ol»oJl  • 

-  Q 

1 

ZmOfL 


1-1 


5  4  C  ,  , 

ol  cl  .  _Ol  "1 

CO  'ojl  "  CO  'ojl 


} 


z-Of  L 


(III-9) 


^  C  c 

Im  i  .y-fs  .a*  -  -  .  -  i  ,a5-c  +  H  .ad..  +  -Si  7*  +  Sk  7*  1 

1  ol  ojl  ol  ojl  ol  15  ol  15  CO  'ojl  CO  OjlJ 

1-1  « 


0  . 


tmOpls 


(III -10) 


Eiqpresslona  for  tbe  coefficients  and  d^^  are  obtained  by  setting  n  ■  0  In 
the  expressions  which  are  given  In  Table  I  of  this  report  [See  Pg.  13  ]. 
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For  th*  lllu*tr«tlv«  problem  under  eonelderetlon,  a  eynmetrleal  fietltioue 
presiure  loading  ebout  the  plane  t  ■  ^  of  the  cylinder  It  applied.  Coneequentlyi 
for  M  ■  Id  bendif  the  ayetea  of  6h  equation*  reduces  to  one  of  31*  ■  3^  equations 
on  the  real  and  Imaginary  portions  of  the  source  strengths. 


A  convenient  form  vhlch  Illustrates  the  order  of  the  matrix  coefficient* 
for  the  source  strength  equations^  vhlch  vas  used  In  the  present  problem  Is  given 
In  Table  V.  The  rows  of  the  matrix,  denoted  by  p,  represent  the  expressions  for 
the  real  cuid  Imaginary  portions  of  the  stresses  and  velocities  at  each  of  the  In¬ 
dicated  field  points  J,  while  the  columns,  denoted  by  q,  indicate  the  coefficients 
of  the  appropriate  source  strengths  at  each  of  the  Indicated  source  points  1.  nie 
last  column  represents  the  forcing  function  vhlch  Is  applied  to  the  cylinder.  It 
should  be  noted  that  the  symbols  d^  which  appear  In  Table  VI  are  defined  In  terms 
of  the  ^  coefficients  of  liable  I  as  follow: 


2H  . 


•^ojlk  *oJ(19 


-l)k 


(Ill-U) 


While  It  Is  not  practical  to  list  all  of  the  numerical  values  of  the  matrix 
coefficients  for  the  illustrative  problem  being  considered,  some  representative 
values  of  these  quantities  for  the  cases  q  ■  1  to  4  and  p  -  1  to  20  are  given  In 
Table  VI  for  the  convenience  of  those  readers  who  may  wish  to  work  through  a 
sanple  cooputatlon. 
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c)  Bcitttlon  cf  th«  S«t  of  SimultaneouB  Llneior  Eguntloin  on  the  Source  Str^ngtht 


01 


*  *01  ®ol 


Th«  solution  of  ths  systsn  of  ths  llnsar  slmultansous  squatlons  on  the 
real  and  imaginary  parts  of  tbs  source  strength  coefficients  can  he  ai>proached 
in  tvo  vays.  Firsts  the  homogeneous  eq[uatlons  In  the  system  can  he  used,  to 
sllalnate  a  nuodber  of  the  unknowns  such  that  the  6n  system  reduces  to  a  series 
of  2N  systems  vhlch  may  then  he  solved  on  an  electronic  computer.  This  procedure 
would  he  Justified  In  a  large  problem  In  %dileh  a  system  of  several  hundred 
equations  would  have  to  he  solved.  For  the  present  problem,  however,  It  Is  a 
comparatively  single  matter  to  solve  the  system  of  equations  directly,  using 
an  IBM  7^90  computer.  The  values  of  the  source  strength  coefficients  are  given 
In  Table  VIZ. 
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Table  VI  -  Numorlcnl  Vnlueo  of  M'ltrlx  CojfflclontB  q  ■  1  to  p  ■ 


1 

2 

3 

4 

p 

1 

-  19.108 

♦  10.552 

♦  3.2568 

♦  1.6250 

2 

+  3.5100 

♦  22.599 

♦  7.7992 

♦  2.4019 

3 

♦.65136 

♦  4.6795 

-  23.672 

+  7.0535 

U 

♦.23214 

♦  1.0294 

♦  5.00B6 

-  20.832 

5 

♦.11003 

♦.41256 

♦  1.1664 

♦  5.2229 

6 

♦.061533 

♦.21305 

♦.48391 

♦  1.24o6 

7 

-.068969 

-.22186 

-.42763 

-.74187 

8 

-.032897 

-.091414 

-.U94i 

-,063712 

9 

-.0085599 

-.016324 

♦.0035093 

♦.070570 

10 

-.0098308 

-.029421 

-.048797 

-.067817 

11 

-.0098062 

-.029349 

-.048678 

-.067653 

12 

-.0097593 

-.029207 

-.048443 

-.067326 

13 

-.0096882 

-.028994 

-.048090 

-.066837 

14 

-.0095939 

-.028712 

-.047623 

-.066189 

15 

-.0095101 

-.028362 

-.047043 

-.065384 

16 

-.0097595 

-.029190 

-.048358 

-.067078 

17 

-.0099396 

-.029730 

-.049253 

-.068334 

10 

-,010031 

-.030002 

-.049705 

-.068963 

19 

-.0016023 

-.0045155 

-.0066488 

-.0077498 

20 

-.0046601 

-.013163 

-.020737 

-.026749 

1  to  20. 


3' 


Tabl*  VZX  •  Source  Strength  Coeffielente 


Source  Point 
1 


ol 

o 

*0. 

(lo®) 

(103) 

(io5)  5., 

-  17. Uu 

37.887 

-  1.7202 

+  .16903 

4  .0040134 

4  .00011319 

-  16.7U1 

31.727 

-  1.3930 

•f  .13424 

-  .043571 

4  .0012951 

•  14.030 

24.660 

-  1.2899 

4  .U856 

4  .12250 

4  .0039739 

.  9.4095 

38.316 

•  1.3449 

♦  .11274 

-  .27677 

4  .0089387 

.  2.9858 

17.601 

-  1.2055 

4  .082023 

4  .43660 

4  .019710 

8.2628 

36.282 

-  .74770 

4  .028792 

•  .81265 

4  .048730 

•  45.238 

7.0447 

•  .65266 

4  .26290 

-  .93899 

•  -13509 

.  47.408 

24.673 

-  1.0528 

♦  .39733 

-  .78204 

-  ..051606 

•  50.463 

28.362 

•  1.2582 

4  .43650 

-  .088151 

-  .025014 
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d)  Evaluation  of  the  Fluid  PreoBure  Field. 


Th«  conqpltx  pressure  and  velocity  con^onents  at  field  points  In  the  fluid 
are  evaluated  from  the  relations 
R 
P, 


>.3 -1 [•  ('ol'.Jl  ♦  Wojl)  *  ‘  (W.J1  -  'ol'ojl)]  <”'■“> 


1-1 

N 


(^ol^ojl  ’  Vojl)  ■  ^  (°0l^0jl  ^ol^Jl)] 


1-1 

R 


[•  (v^ji  -  5oi“4i)  -  ^  ^  "^01^0, i)] 

1-1 


(111-13) 


(III-lJv) 


The  coefficients  and  their  space  derivatives  are  evaluated  from 

Eqs.  (A-1)-(A-38)  of  ^pendlx  A  of  this  paper  vhere  R  Is  the  distance  between 
the  source  point  1  and  the  field  point  J  In  the  acoustic  fluid  at  which  the 
pressure  and/or  velocities  are  being  coniputed> 

£-a  "o  -  ®o  - 

The  quantity  X  -  •-*-  Is  given  by  Eq,  (A-lt)  where  p  ■  cos  t  i  I  ■  sin  C/  snd 
R^  represents  the  distance  from  the  field  point  J  at  which  the  pressure  is  being 
evaluated,  to  a  reference  point  0  which  for  convenience  Is  located  at  the  center 
of  the  cylinder,  l.e.  at  the  point  0  defined  by  the  coordinates  r  ■  0  and  s  -  ^  , 
yig  (II1-3). 


The  far  field  pressures  In  the  fluid,  l<e«  the  pressures  at  distances  R 

irtilch  are  large  multiples  of  the  cylinder  dimensions,  are  of  practical  interest. 

A  sliiq>llfled  asyoiptotlc  expression  for  the  fluid  pressure,  p^  .  In  the  far  field, 

o,j 

was  derived  In  Reference  [1]  In  terms  of  the  fluid  source  strength  coefficients 
and  a  prescribed  distance  and  slope  angle,  R^  and  C  respectively,  as  shown 
In  Fig.  (ZZI-3)<  P'or  the  axl symmetrical  loading  and  geometry  of  the  present 


problemi  the  asymptotic  value  of  the  far  field  fluid  pressure  Is  given  by  the 


relation 


Ico  R. 


'oj 


r  -  ^*'3*  T  /V  -  lo’ti 

•  “2“  *  *^0  ^3*"  ****  "  *  J 

1-1 


(III-15) 


It  should  be  noted  that  once  the  pressure  p  ,  has  been  evaluated  from  Eq«  (III>15} 

for  a  given  point  P.  on  the  ray  specified  by  R  and  (  ,  it  con  easily  be  evaluated 
J  ® 

for  any  other  point  on  the  ray  by  changing  the  scale  factor 

lo)  R^ 

0 


since  the  quantity  In  the  bracket  will  be  a  constant  for  a  given  ray.  nie  use 
of  Eq.  (XIZ«15)  in  the  evaluation  of  the  far  field  fluid  pressure  greatly 
slnqillfles  the  confutation  vhlcb  vould  be  required  If  Eq.  (III>12)  vere  to  be 
usedi  since  it  eliminates  the  necessity  for  the  numerical  computation  of  the 
^ojl  coefficients  for  each  field  point  J  In  the  fluid. 


A  slsfle  formula  for  the  evaluation  of  the  absolute  value  of  the  pressure 
p.  ,  in  the  far  field  la  given  by  the  relation 


-  -5" 


IV  ik^s  sin  C 

Z®oi'  K -  V^3?cost) 


i-1 


®c  .  ^ 

%  • 


ry 

LrJ 


(111-16) 


where  refers  to  sane  reference  point  on  the  rey  under  consideration  and  R^ 
refers  to  the  field  point  on  the  ray  at  which  the  pressure  Is  to  be  evaluated. 
For  points  which  are  closer  to  the  cylinder  so  that  the  asynf totlc  formula  for 
the  fluid  pressure  cannot  be  used,  the  absolute  value  of  the  pressiire  .  is 
given  by  the  relation, 


(111-17) 


l-I  i-1 


Ojl 


Vojl 


Rum«rleal  computations  for  the  abaoluto  value  of  the  pressure  p 


o»J 


have 


been  carried  out  along  rays  ranging  from  0®(224®)90®  as  shown  in  Figs.  (XII-U) 
and  (III-S).  In  each  figure,  the  lines  of  constant  pressure  in  the  fluid  are 
shown.  The  results  which  ore  plotted  have  all  been  derived  using  Eq.  (1II«17) 
and  checked,  beyond  -  20n,  by  Eq.  (III-I6).  For  the  benefit  of  those  who  may 
wish  to  work  through  a  sample  confutation,  a  representative  set  of  numerical 
values  is  given  in  Table  VIXI.  The  values  of  the  constants  appearing  in 
Eq.  (XXX>l6)  are  given  in  Table  XX.  Xn  the  present  example,  Eq.  (XXX>l6)  can 
be  used  with  good  accuracy  for  those  points  for  which  R^  >  20a. 


It  must  once  again  be  enphasized  that  these  numerical  results  may  not 
represent  an  accurate  solution  of  this  problem  because  too  few  bonds  were 
purposely  used  in  the  computational  work  for  the  cylinder  and  fluid  source 
strength  coefficients.  This  was  done  In  order  to  present  the  simplest  possible 
numerical  example  which  lUustrated  all  the  conplexitles  of  the  application  of 
the  general  method.  For  a  practical  problem  with  a  real  transducer,  a  considerably 
finer  finite  difference  breakup  would  be  required  in  order  to  evaluate  accurate 
values  of  the  pressures  in  the  fluid. 


I 

I 
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Tabl«  VIII  »  Computation  for  Prescure  In  Fluid  -  Numerical  Extunplei 

R 

Roy  3#  i  -  »»5®,  -  7.07. 


■i 

— 

• 

* 

■ 

■ 

B 

t 

’'oji 

■ 

1 

1 

•  .006162 

♦  .0019395 

B 

i 

1 

-  .017912 

+  .0056330 

1 

1 

-  .027798 

+  .0087756 

fl 

1 

-  .034809 

*  .011021 

5 

i 

1 

-  .036191 

.012058 

6 

u 

I2 

1 

-  .037537 

.011626 

■ 

1 

1 

-  .032898 

♦  .018512 

B 

1 

-  .022013 

♦  .031663 

B 

1 

1 

•  .0054818 

*  .038822 

10 

1 

1 
’  5 

♦  .013308 

+  .037387 

11 

1 

•5 

+  .029733 

+  .026771 

12 

1 

-1 

+  .03918 

+  .0089625 

13 

u 

15 

-  1 

♦  .040008 

-  .0011364 

14 

1 

•  1 

♦  .037455 

-  .00062594 

15 

-  1 

♦  .032249 

-  .00022219 

16 

•  1 

*  .024762 

*  .000022490 

17 

1 

Ti 

-  1 

♦  .015590 

*  .000097948 

18 

1 

15 

-  1 

*  .0053192 

+  .000047970 

^o,J  ■  •  'OOlMDSl  j  -  .002333  pii. 


ts  [R  ■  20a]. 

c _  _ 


^ay  tColculated  from 
Eq.  (III-17)  at  R  -  20al 

0.910 

0.883 

0.824 

0.788 


0.798 


IV  DiGciTSSion  ACT)  concLusr.ious 


Thle  paper  la  port  of  a  study  on  tha  developraent  of  methodo  for  treating 
the  forced  rlbratlona  of  elastic  bodies  of  revolution  oubinorccd  In  an  acoustic 
fluid.  A  method^  based  on  a  potential  theory  approach,  for  the  evaluation  of 
the  pressure  field  In  an  Infinite  acoustic  fluid  due  to  the  harmonic  excitation 
of  an  elastic  circular  cylindrical  body  of  finite  length  has  been  derived  In 
Reference  [1]  and  Is  illustrated  In  the  present  paper. 

The  Illustrative  example  Is  presented  for  the  case  of  an  axlsynmetrleal 
(•  independent)  excitation  of  the  cylinder.  The  application  of  the  method  to 
more  general  excitations  of  the  cylindrical  body,  vhlch  vary  In  both  s  and  9, 
follow  a  similar  conputatlonal  pattern  and  are  described  In  Port  XI  of  this 
report. 

The  general  procedures  nay  find  application  In  problems  relating  to  the 
response  of  fluid  to  the  harmonic  excitation  of  large  cylindrical  transducers. 
In  addition,  the  general  approach  Is  also  being  applied  to  problems  Involving 
the  evaluation  of  fluid  pressures  produced  by  the  harmonic  excitations  of 
submerged  cylindrical  shells  of  finite  length. 


FI6.M-3  GEOMETRY  FOR  THE  EVALUATION 
OF  FLUID  PRESSURES 


FIG  m-4  FLUID  PRESSURES  i  Ro/i 
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FIG. 1-5  FLUID  PRESSURES  J  Ro/ 
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APPENDIX  A  -  EVALUATION  OF  THE  REAL  AND  IMAGIKARY  PORTIOIIS  OF  TIIE  COEFFICITCirrS 
"•njl  *  ^njl  *  DERIVATIVES. 


Por  computational  purposea,  the  complex  coefficients 
are  written  In  terns  of  their  real  and  Imaginary  ports,  e.g. 


“njl  -  »njl  ♦  ^  “njl 


(A-1) 


The  following  nomenclature  which  Is  similar  to  that  used  In  Appendix  B 
of  Reference  [1]  for  the  case  n  ■  0,  Is  used  In  this  Appendix: 


•  2  2  2  • 

R  •  ((s-s)  ♦  r  ♦  r  -  2ir  cos  t) 


i 


(A-2) 


.  25  ,  1  ■  1,2,3  (a.3) 

S  Cj 

X  -  k^R  -  [(e-D®  ♦  5®  +  p®  -  2pp  cos  yji  (A-U) 

where 

m  m 

l-|>  l"|i  P"J  P"|  (A-5) 


Specific  equations  for  the  coefficients  and  their  space 

derlwatlves  will  be  given  In  the  Appendix.  The  corresponding  ejqsresslons 
for  the  coefficients  and  7^^  and  their  derivatives  are  obtained  by 
replacing  ^2  "  ^njl^  <>*“  ^3  ■  T  ’'njl^  respectively 

in  the  equations  which  follow. 


tvrt  of  Coefficient  >nd  PeriratlYe* 


•)  IjLi 


®.0t  •  ■  ^  »t  M 


“Jjl  “  J  (f-»  “•  CM  »♦  4t 


%«‘5ri  (l-l)»i(*)  CM  nt  «t 
0 


•■^1  "  ^/  *"•  t )%(♦)!  «0»  >»♦  4t 

0 

*^i  •  ^ 

0 

J 

^ /t3(e-c)»2(*)*l«-i>^3<*>i  «t 

C 


(A-6) 


(A-7) 

{A-8) 

(A-9) 

(A-10) 

(A-11) 

(A-12) 


fi?P  ^ 


«*•  PO0  nt  At  (A.13) 


l(p-0  eol  t)tg(<')*{l-C)®(p-P  co»  ♦)>2(t)l  oos  nt  df  (A-lU) 


f  (t)  .  ^.1  X 

(A-15) 

f  (t)  .  1] 

(A-16) 

f  _  j{15-ex*)co«  X  ♦  (15-X®)  X  sin  Xj 

3  vT 

<A-17) 

M  kjLA 


For  the  ease  i  ■  l«e>  vhere  the  field  point*  and  the  source  points  lie 
on  the  same  bend^  special  formulas  are  required  for  the  computation  of  , 

^nll  '  ^nll  their  space  derivatives.  The  reader  is  referred  to  Appendix  C 
of  Reference  [1]  for  the  derivation  of  these  expressions  for  the  case  n  ■.  0. 

The  corresponding  formulas  for  n  ^  0  are  given  belov. 

c-  r  w  r  sin  (~)n 

“nil  •  ’  kj[  ^  [ - 

nt  di|r  J  on  surfaces  z«0^L 

(A-19) 

nt  dt  +  ^  J  on  surface  r-a 


•^Sll  " 


3 

i. /‘o  -p  cos  ♦)  cos 


-P  cos  t) 


cos 


0 


J  on  surfaces  e-O,! 
J,  on  surface  i>a 


(A-20) 


-2«rr 
•  ®nll 


i  on  surface  r^a  (A-21} 
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«  COB  (~yn 

nt  dt  ♦  [-  j  ®“ 


A- . « 


>^L‘i  ■ « 


i-«rr  rt 


9 


-  It  J  ^8^*^  ®®*  *'•' 


J  on  turf aces  taO,L 

(A-23) 

J  on  all  surfaeea 

(A-2U) 

J  on  all  surfaces 

(A-25) 

At  J  on  all  surfaces 

(A-26) 
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Imaginary  Port  of  Coefficient  and  Derlvativee. 


c)  For  all  combinations  of  1  and  J. 


r 

2"  J 


Bln  X 


*njl  ■  2k  J  X 
0 


COB  nt  dt 


iA-27) 


k^p  * 


aa^l  *2^/  ®°® 


(A-28) 


* 


“«i  ■  /«•*>  \<*>  “•  “♦ 


(A-29) 


it?5  * 


(A-30) 


k^R  * 
*iP 


^  J  *  U-i)^ 


(A-31) 


k?P  " 


,3a”*  -  ~  J  [3{l’l)  FgCt)  +  (l-l)^  FjCt))  COB  ni  dt 


(A-32) 


2«rt 

•®njl 


R^p  " 


gi- J^(P-P  co»  FgCt)  COB  Ulf 


dy 


(A-33) 


ii?P  * 

^ (t*?)(p-P  COB  t)^  Pjdr)]  c^.g  ni|f  df 


(A-34) 
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.3>rxc 
la 


J l(p-2  eo»  ♦)  ^2^1)  ♦  (t-l)®(p-p  POP  ♦)  POP  nt  dp  (A.35) 


vhere 


‘3+20'B1Io  ••• 

Fx  cos  X  -  Bln  xl 

L - ;s — J  - 


X  <  0.06 

X  >  0.06 


(A-36) 


•  1a 

FgU)-  15-^0^7555  *** 

{3-X^)  sin  X  -  3X  cos  X 

^  ,  "■  .  top 


X  <  0.40 


X  >  o.4o 


(A-37) 


snd 


1  .  X®  X" 


*'3(7)-  -  io5  *  IS9  ■  BSaSS  — 

(6X^«  15)  Bln  X  (15  -  X^)  X  cos  X 


X  <  1.0 


X.>  1.0 


(A-30) 


■  *  a 

The  quantities  and  and  their  derivatives  do  not  contain  an 

infinite  discontinuity  in  their  integrands  at  the  value  defined  by  R  ••  0. 
Consequently^  Eqs.  (A«27}'>(A>35)  nBy  be  evaluated  directly  and  there  is  no 
need  for  special  formulas  as  in  the  case  of  the  real  (single  barred)  components. 


APPEWDIX  B  -  DIAGRAMMATIC  SUMMARY  OF  COMPUTATIONS. 

MCtiBAry  eonputatlons  for  toe  •valuation  of  tha  prasaur*  fiald  la 
tba  fluid  dua  to  tha  tina  hamonle  axoltatloa  of  an  alaatlc  cylindrical 
body  of  finite  length  ara  suanMuriaad  in  Tigs.  (B-l)  to  (B-3)*  Tha  conq^utationa 
ara  convaniantly  hrokan  into  thraa  baaic  blocks  aa  ahown  in  tba  diagrams. 

Tha  order  of  cooputatlons  ean>  if  nacassary,  ba  raorganisad  in  accordance  vith 
tha  specific  raquiraaents  of  the  particular  electronic  con^utar  being  used. 
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FIO.  B-2  -  DIAGRAMMATIC  SUMMARY  -  BLOCK  II. 


Computation  for  unknown  source 
Strengths  C^^  , 
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no.  B-3  -  DIAGKA1..L\TIC  SUIilAl^Y  -  BLOCK  III. 
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Code  ••21 
Code  ••23 
Code  ••aj 
Code 

Code  ••••a 
Code  ••43 

Code  689,  Mr.  I.  Cook 
Code  1500 

Chief,  Bureau  of  Navol  Weapons 
Deportment  of  the  Navy 
Washlncton  25,  D.C. 

Attn:  DAAD,  Alrfrome  Dcolcn 
DLI-3,  Tech.  Library 
R-12,  Clilef  Scientist 
RMOA,  IX}  &  Alrfrtuno  Br. 

RU,  ASW  Division 
RRRE,  Research  Br, 


Cpcclf'l.  Iroji.'ctii  Office 
Durcf  u  of  n.:.v  1.  './  •';»ou.'j 
Dopartniont  of  tlie  H.ivy 
Wuuhln,'',ton  25,  D.C. 

Attn:  Code  SI’-OOl,  Chief  Scientist 
Code  sr-20,  Tech.  Director 

Chief,  Durociu  of  Yards  &  Docks 
D  jm'tiii'ut  of  the  Navy 
W:ird)lr',;l-.on  25,  D.C. 

Attn:  Co>1'!  JO,  P.eBcarch 

Code  Tech.  Ubraxy 

Commindlns  Officer  &  Director 
David  Tay).or  Model  Benin 
Waohlnf'ton  7,  D.C. 

Attn:  Code  lOO,  Mr,  R.T,  McGoldrlek 
Code  lOoD,  Dr.  M.  Straasberg 
Code  l40,  Tech.  Info.  Dlv. 

Code  53O,  Dr.  F.  Thcilhelmer 
Code  5^3,  ^tr.  A.O.  Sykes 
Code  700,  Dr.  A.H.  Kell 
Code  720,  Mr,  E.E.  Johnson 
Code  731,  Mr,  J.C.  Pulos 
Code  740,  Dr.  W.J.  Bette 
Code  760,  I'tr.  E.  Noonan 
Code  761,  Dr.  E.  Buchman 
Code  771,  Dr.  R^  UobowltK 

DTI®  Undenmter  Explosion  Res,  Dlv, 
Norfolk  Naval  Shipyard 
Portemouth,  Vlr(.;lnlB 
Attn;  Mr.  D.S.  Coben 

Commander 

U.3.  Naval  Ordnance  Laboratory 
White  Oak,  Maryland 
Attn:  HL,  Technical  Library 
D,  Technical  Director 
OU,  Underwater  Weapons 
RA,  Acouatlco  Division 

Commanding  Officer  &  Director 

U.8.  Navy  Underwater  Sound  Laboratory 

Port  Trumbull 

New  London,  Connecticut 


Cowiivmdln'j  Officer 

U.Di  Itwil  Kdno  Dofonoa  Laboratory 

rciKJiia  City,  Florida 

CojiiiiKtiidor 

U.S.  Nnvnl  Air  Dovclopwent  Center 
Johnavllle,  I’ennoylvunlo 

Director 

U.S.  Navy  Underwater  Sound 
Reference  Laboratory 
Office  of  Naval  Rconcirch 
P.O.  Box  8337 
Orlando,  Florida 

Commanding  Officer  &  Director 
U.S.  Navy  Electronics  Laboratory 
Sen  Diego  52,  California 

Commander 

Portsmouth  Noval  Shipyard 
Portsmouth,  New  Hampshire 

Commander 

Mure  Island  Naval  Shipyard 
Vallejo,  California 

Director,  Materials  Laboratory 
Rev  York  Naval  Shipyard 
Brooklyn  1,  New  York 

Offleer«ln-Chargs 
Naval  Civil  Engineering  Research 
&  Evaluation  Laboratory 
U.S.  Naval  Construction 
Battalion  Center 
Port  Hueneme,  California 

Director 

Naval  Air  E]q|)erlment  Station 
Naval  Air  Material  Center 
Naval  Base 

Philadelphia  12,  Pennsylvania 
Attn I  Structures  Laboratory 

Of f 1 cer-ln-Charge 

David  Taylor  Model  Basin 

Underwater  E]q)loslon  Research  Division 

Norfolk  Naval  Shipyard 

Portsmouth,  Virginia 

Attn I  Dr.  H.N.  Schauer 


Conuiinndcr 

U.S.  Naval  Proving  Ground 
Dahlgrcn,  Virginia 

Bupcrvlcor  of  Shipbuilding,  U6N 
and  Navi  I  liu.poctor  of  Ordnance 
Genortil  Dynur.iicn  Corporation 
Electric  IJoat  Division 
Oroton,  Connecticut 

Comma  ndur 

Naval  Ordnance  Test  Station 
China  Lake,  California 
Attni  Physico  Division 
Mechanics  Division 

Commanding  Officer 
Naval  Ordnance  Teat  Station 
Underwater  Ordnance  Division 
3202  E.  Foothill  Boulevard 
Pasadena  8,  California 

Commanding  Officer  &  Director 
U.S.  Naval  Engineering  E}q>eri> 
ment  Station 
Annspolls,  Maryland 

Superintendent 

U.S.  Naval  Postgraduate  School 
Monterey,  California 

F,  Air  Force 

Commander 

Air  Material  Command 
Wrlght-Patterson  Air  Force  Base 
Dayton,  Ohio 
Attni  MCREX-B 

Structures  Division 

Commander,  WADD 

Vrlght'Patterson  Air  Force  Basa 
Ohio 

Attni  WWRC 

WWRMD8 

WHBKDD 

Director  of  Intelligence 
Headquarters,  U.S.  Air  Force 
Washington  25,  D.C. 

Attni  P.V.  Branch  (Air  Targets  Dir) 


Coii'jl'M’.iidcr 

Air  Korco  Office  of  Scientific 
Re soar ch 

VELohlneton  25,  D.C. 

Attni  Mochanlce  Dlvlolon 

0*  Other  Oovernmont  Activities 

U«S>  Atomic  Energy  Cotnmleelon 
Woohlngton  25,  D.C. 

Attn:  Director  of  Research 

Director 

National  Bureau  of  Standard! 
Waohlngton  25,  D.C. 

Attn:  Division  of  Mechanics 

National  Aeronautics  &  Space  Adn. 
1512  H  Street,  N.W. 

Washington  25,  D.C. 

Attn:  Chief,  Dlv.  of  Research 
Information 

Director 

National  Aeronautics  A  Space  Adm. 
Langley  Research  Center 
Langley  Field,  Virginia 
Attn:  Structures  Division 

H.  Contractors  and  Others 

Dr.  M.Li  Baron 

Paul  Weldllnger,  Consulting  Engineer 
770  Lexington  Avenue 
New  York  21,  New  York 

Professor  H.H.  Blelch 
Department  of  Civil  Engineering 
Columbia  University 
6l0  Mudd  Building 
New  York  27,  Now  York 

Profeesor  B.A.  Boley 
Department  of  Civil  Engineering 
Columbia  University 
610  Mudd  Building 
New  York  27,  New  York 


Dr.  M.A.  Brull 

Englnuurlng  lUch.inlca  Division 
Unlvcrolty  of  fi-nnHylvuiila 
Philadelphia  4,  Pennsylvania 

Dr.  F.L,  DlM.if;;lo 

Dept,  of  Civil  Engineering 

Columbia  University 

610  Kudd  Building 

Now  York  27,  New  York 

Profeesor  D.C.  Drucker 
Division  of  Engineering 
Brown  University 
Providence  12,  Rhode  Island 

Dr.  Ira  Dyer 

Dolt,  Beranek  &  Newman,  Znc. 

50  Moveton  Street 
Cambridge  3^,  Massachusetts 

Professor  A.C.  Erlngen 

Dept,  of  Aeronautical  Engineering 

Purdue  University 

Lafayette,  Indiana 

Profeesor  J.N.  Ooodlcr 
Dupnrtment  of  Mechanical  Engineering 
Stanford  Unlvi'mlty 
Stanford,  California 

Dr.  Joeh  E.  Greonopon 

J.O.  Engineering  Research  Assoc. 

3709  Cnllnwoy  Avenue 
Baltimore  15,  Maryland 

ProfesBor  Philip  0.  Hodge 
Armour  Research  Foundation 
10  West  35th  Street 
Chicago,  Illinois 

Professor  Nicholas  J.  Hoff 
Dept,  of  Aeronautical  Engineering 
Stanford  University 
Stanford,  California 


Dr.  C.W,  Jlorton 
Defence  Rooeorch  L:iboratory 
Unlvorolty  of  Texas 
Austin  12,  Texas 

Dr.  M.C.  Juncer 

CcunbrldQQ  Acoustical  Associates  ' 
129  Mount  Auburn  Street 
Ceunbrldce  ^taosachueetts 

Professor  Joseph  Kcn^ner 
Dept,  of  Aeronautical  Engineering 
&  Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
333  Jay  Street 
Brooklyn  1,  New  York 

Professor  J.M.  KLosner 
Dept,  of  Aeronautical  Engineering 
&  Applied  Kechonles 
Polytechnic  Institute  of  Brooklyn 
333  Jay  Street 
Brooklyn  1,  New  York 

Professor  E.H.  Lee 

Brown  University 

Division  of  Applied  Mathematics 

Providence  12,  Rhode  Island 

Professor  R.D.  Mlndlln 
Dept,  of  Civil  Engineering 
Columbia  University 
618  ^kldd  Building 
New  York  27,  New  York 

Professor  P.M.  Noghdi 
University  of  California 
College  of  Engineering 
Berkeley  4,  California 

Professor  N.M.  Newnork,  Head 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbona,  Illinois 

Professor  F.  Pohle 
Department  of  Mathematics 
Adelphl  College 
Carden  City,  New  York 


Professor  W.  Pri’"?r,  Chairman 
Phyolcrl  Scloncjc  Council 
Brown  University 
Providence  12,  Rhode  Island 

Professor  F.V.  Romano 
Dept,  of  Aeronautical  Engineering 
&  Applied  M-chnnlcs 
Polytechnic  Institute  of  Brooklyn 
333  Joy  Street 
Brooklyn  1,  New  York 

Professor  E.  Reiss 

Institute  of  t'-ithcmatlcal  Sciences 

New  York  University 

25  Waverly  Place 

Hew  York  3»  Ucw  York 

Dr.  V.L.  Salerno 

Don  Dosco  Institute  for  Research 
Rumsej',  New  Jersey. 

Professor  A.S.  Veletsos 
Deportment  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois 

Brown  University 
Research  Analysis  Croup 
Providence,  Rhode  Islond 

Electro  Nuclear  Systems  Corporation 
Technical  Reoenreh  Division 
6001  Norfolk  Avenue 
Bethesda  l4,  Maryland 

Hudson  Laboratories 
Columbia  University 
l4^  Palisades  Street 
Dobbs  Ferry,  Now  York 

Lament  Ccologlcal  Observatory 
Columbia  University 
Torre  Cliffs 
Palisades,  New  York 

Ordnance  Research  Laboratory 
University  of  Pennsylvania 
University  Park,  Pennsylvania 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  Massachusetts 
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